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Abstract 

We report on the evidence for the decays B°—>-DfiT~ and B°^Dj K + and the results of a search 
for the decays B°^D* + tt~ and B°^D*~ K + from a sample of 84 million T(AS) decays into B 
meson pairs collected with the BABAR detector at the PEP II asymmetric-energy e + e~ collider. The 
measured B°^-Dfir~ yield has a probability of less than 10~ 3 to be a fluctuation of the background 
and we measure the branching fraction B{B°^D+^) = (3.2 ± 0.9 (stat.)rbl.O (syst.)) xlO" 5 . The 
measured B°^D~ K + yield has a probability of less than 5 x 10~ 4 to be a fluctuation of the 
background and we measure the branching fraction B(B°^D~ K + ) = (3.2 ± 1.0 (stat.)il.O (syst.)) 
xl0~ 5 . We also set 90% C.L. limits B(B°^D* + ir-)< 4.1xl0~ 5 and B(B°^D*-K + )< 2.5xl0~ 5 . 
All results are preliminary. 
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The measurement of the CP-violating phase of the Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix ||] is an important part of the present scientific program in particle physics. CP violation 
manifests itself as a non-zero area of the unitarity triangle Q. While it is sufficient to measure 
one of the angles to demonstrate the existence of CP violation, the unitarity triangle needs to be 
overconstrained by experimental measurements, in order to demonstrate that the CKM mechanism 
is the correct explanation of this phenomenon. Several theoretically clean measurements of the 
angle (3 exist ||, but there is no such measurement of the two other angles a and 7. A theoret- 
ically clean measurement of sin(2/3 + 7) can be obtained from the study of the time evolution of 
the B°— >.t)(*)~7r + Q decays, of which a large sample is already available at the B-factories, and 
of the corresponding Cabibbo suppressed mode B°^D^ + tc^ This measurement requires the 
knowledge of the ratio between the decay amplitudes = \A(B ^D^ + Tr-)/A(B ^D^-TT + )\. 
Unfortunately the measurement of \A(B°— >D^ + tt^)\ via the measurement of B(B°— >D(*^ + ir~) is 
not possible with the currently available data sample due to the presence of the copious background 
from B°— >Z)(*) + 7r~. However, we can measure B(B°— >Di*^ + ir~) and relate it to using SU(3) 

symmetry: R^ 2 oc ^^ "^^ <i) _ 7r + | , where the proportionality constant is, to first approximation, 

the ratio of the D$ and the decay constants (5|. The decays B°— >Z)i*' )+ 7r~ have also been 

proposed to be used for the measurement of | V u b/V c b\ ||. 

The decays B°^D { S * ] ~K + are a probe of the dynamics in B decays because they are expected 
to proceed mainly via a W-exchange diagram, not observed so far. In addition, theses modes can 
be used to investigate the role of final state rescattering since its presence can substantially increase 
the expected rates [ffl. In this letter we present measurements of the branching fractions for the 
decays B°^D { S * )+ -K~and B°^D { s * y K + . 

This analysis uses a sample of 84 million T(4S') decays into BB pairs collected in the years 1999- 
2002 with the BABAR detector at the PEP-II asymmetric-energy B-factory f|. Since the BABAR 
detector is described in detail elsewhere |J, only the components of the detector crucial to this 
analysis are summarized below. Charged particle tracking is provided by a five-layer silicon vertex 
tracker (SVT) and a 40-layer drift chamber (DCH). For charged-particle identification, ionization 
energy loss (dE/dx) in the DCH and SVT, and Cherenkov radiation detected in a ring-imaging 
device are used. Photons are identified and measured using the electromagnetic calorimeter, which 
comprises 6580 thallium-doped Csl crystals. These systems are mounted inside a 1.5 T solenoidal 
superconducting magnet. We use the GEANT [1C] software to simulate interactions of particles 
traversing the BABAR detector, taking into account the varying detector conditions and beam 
backgrounds. 

We select events with a minimum of three reconstructed charged tracks and a total measured 
energy greater than 4.5 GeV as determined using all charged tracks and neutral clusters with energy 
above 30 MeV. In order to reject continuum background, the ratio of the second and zeroth order 
Fox- Wolfram moments [jll| must be less than 0.5. 

So far, only upper limits on the modes studied in this letter exist [JDJ. Therefore the selection 
criteria are optimized to maximize the ratio of signal efficiency over the square-root of the expected 
number of background events. 

The Df mesons are reconstructed in the modes Df—xpir^, K®K + and K*°K + , with cj)^K + K~~ , 
- and K*°^K-ir + . The K° s candidates are reconstructed from two oppositely-charged 
tracks with an invariant mass 493 < M % + % - < 501 MeV/c 2 . All other tracks are required to originate 
from a vertex consistent with the e + e~ interaction point. In order to identify charged kaons, two 
selections are used: a pion veto with an efficiency of 95% for kaons and a 20% pion misidentification, 
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and a tight kaon selection with an efficiency of 85% and 5% pion misidentification probability. If 
not otherwise specified, the pion veto is always adopted. The 4> candidates are reconstructed from 
two oppositely-charged kaons with an invariant mass 1009 < M K + K - < 1029MeV/c 2 . The K*° 
candidates are constructed from the K~ and a ir + candidates and are required to have an invariant 
mass in the range 856 < M K - W + < 936MeV/c 2 . The polarizations of the K*° 0) mesons in the 
Df decays are also utilized to reject backgrounds through the use of the helicity angle 9h, defined 
as the angle between one of the decay products of the K*° (cj)) and the direction of flight of the 
meson itself, in the meson rest frame. Background events are distributed uniformly in cos 9h 
since they originate from random combinations, while signal events are distributed as cos 2 Oh- 
The K*° candidates are therefore required to have |cos0#| > 0.4, while for the <f> candidates we 
require | cos##| > 0.5. In order to reject background from D + — >K® tt + or K*° 7r + , the K + in the 
reconstruction of Df^K® K + or K*° K + is required to pass the tight kaon identification criteria 
introduced above. Finally, the Df candidates are required to have an invariant mass within 10 
MeV/c 2 of the nominal mass |12| |. 

We reconstruct D* + candidates in the mode D* + '— >L>+7, by combining and photon candi- 
dates. Photons that form a 7r° candidate, with 122 < M 77 < 147MeV/c 2 , in combination with any 
other photon with energy greater than 70 MeV are rejected. The mass difference between the D* + 
and the Df candidate is required to be within 14 MeV/c 2 of the nominal value fl^| . 

We combine Df or D* + candidates with a track of opposite charge to form B°— >D^ K + 

r\ I 

or B —>Ds vr _ candidates depending on whether they pass the tight kaon selection criteria. In 
order to reject events where the Df comes from a B candidate and the pion or kaon from the other 
B, we require the two candidates to have a probability greater than 0.25% of originating from a 
common vertex. The remaining background is predominantly combinatorial in nature and arises 
from continuum qq production. In order to suppress it using the event topology, we compute the 
angle (9t) between the thrust axis of the B meson decay product candidates and the thrust axis 
of all the other particles in the event. In the center-of-mass frame (cm.), BB pairs are produced 
approximately at rest and produce a uniform cos 6t distribution. In contrast, qq pairs are produced 
back-to-back in the cm. frame, which results in a | cos 9t\ distribution peaking at 1. Depending 
on the background level of each mode, | cos 9t\ is required to be smaller than a value which ranges 
between 0.7 and 0.8. We further suppress backgrounds using a Fisher discriminant T constructed 
from the scalar sum of the cm. momenta of all tracks and photons (excluding the B candidate 
decay products) flowing into 9 concentric cones centered on the thrust axis of the B candidate [ 14 1 . 
The more spherical the event, the lower the value of T . We require T to be smaller than a threshold 
which varies from 0.04 to 0.2 depending on the background level. 

We extract the signal using the kinematic variables toes = J E^ 2 — P* ) 2 an d AE = 



J2i y m f + Pi 2 — E£, where E£ is the beam energy in the cm. frame, p* is the cm. momen- 
tum of daughter particle i of the B meson candidate, and nii is the mass hypothesis for particle i. 
For signal events, tties peaks at the B meson mass with a resolution of about 2.5 MeV/c 2 and AE 
peaks near zero, indicating that the candidate system of particles has total energy consistent with 
the beam energy in the cm. frame. The AE signal band is defined by \AE\ < 36 MeV and within 
it we define as signal candidates the events with ttjes > 5.27GeV/c 2 . 

After the aforementioned selection, three classes of backgrounds remain. First, the amount 
of combinatorial background in the signal region is estimated from the sidebands of the tties dis- 

dN _ „ . /-. _ n .n 7^*2 qvti Lt/1 _ ~2 / J?*1\ 



tribution and is described by a threshold function = xJl — x 2 /E^ 2 exp [— £ (l — x 2 /E£ 2 )] , 



characterized by the shape parameter £ [15]. 
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AE(GeV) 

Figure 1: The AE distribution in data compared with the distribution in the combinatorial back- 
ground, estimated from the m^s sidebands, and with the cross-contamination, which is estimated 
from the M^™ d sidebands. The insert shows separately the AE distribution of the contributions 
to the cross contamination as expected from the simulation. The reflection background is normal- 
ized to the known branching fractions |12| , while the normalization of the charmless background is 
arbitrary. 



Next, B meson decays such as B°— >D + ir~, p~ with D + —>Kgit + or K*°tt + can constitute a 
background for the B°^-Dfir~ mode if the pion in the D decay is misidentified as a kaon (re- 
flection background). This background has the same tties distributions as the signal but different 
distributions of AE. The corresponding background for the B°^>D~ K + mode (B°^D~ K + , K* + ) 
has a branching fraction ten times smaller. Finally, rare B decays into the same final state, such 
MS 7r or K^°K + K (charmless background), have the same mES an d AE distribu- 

tions as the signal. Figure shows the AE distribution for the signal and for the various sources of 
background. 

The branching fraction of the charmless background is not well measured and we therefore need 
to estimate the sum of the reflection and charmless background (referred to as cross- contamination) 
directly on data. This is possible because both of these backgrounds have a flat distribution in 
the D+ candidate (M%™ d ) mass while the signal has a Gaussian distribution. Possible peak- 
ing background from B^DfX decays is negligible, as determined from simulation. The cross- 
contamination to the decays B°^D* + it~ and B°^D*~ K + is dominated by the reflection back- 
ground which we estimate from simulation. Cross-feed between B°— >Di*^ + ir~ and B°— >D^ K + 
modes has been estimated to be less than 1%. 
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5.25 5.3 5.25 5.3 

m ES (GeV/c 2 ) 

Figure 2: The m E s distributions for the B°^D+tt- (top left), B°^D* s + ir~ (top right), B°^D-R + 
(bottom left), and B®^D* S ~ K + (bottom right) candidates on data after the selection, within the 
AE band. The fits used to obtain the signal yield are described in the text. The contribution from 
each mode is shown. 



Figure shows the toes distribution for each of the modes in the AE signal band. We perform 
an unbinned maximum-likelihood fit to each toes distribution with the threshold function to char- 
acterize the combinatorial background and a Gaussian function to describe the sum of the signal 
and cross-contamination contributions. The mean and the width of the Gaussian distribution are 
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Table 1: The number of signal candidates (N s i g t, ox ), the Gaussian yield (N gaus ) and the com- 
binatorial background (N corn i,) as extracted from the likelihood fit, the reconstruction efficiency 
(e), the cross-contamination (N cross ), the probability (Pbckg) of the data being consistent with the 
background fluctuating up to the level of the data in the absence of signal, the measured branching 
fraction (£>), and the 90% confidence level upper limit. N gaus , N com f, and B are not available for 
modes with too few events. N cross is not reported if no event is found in the mass sideband. 



B mode 


N sigbox 


N gaus 


N comb 


cross 


e(%) 


-Pbckg 


B{10- b ) 


90% C.L. 

(io~ 5 ) 


B°^DTtt- 






















9 


8.0 ±3.0 


2.1 ± 0.7 


< 0.7 


16.9 


1.4 x 10" 


3 


3.1 ± 1.2 


- 


Dt->K*°K+ 


12 


9.2 ± 3.4 


3.8 ± 1.0 


2.9 ± 1.8 


9.6 


2.3 x 10" 


2 


3.5 ± 1.9 


- 


D+-*K° S K+ 


5 


4.2 ± 2.2 


1.9 ±0.6 


1.2 ± 1.4 


12.3 


8.3 x 10" 


2 


2.4 ± 1.8 


- 


all 


26 


21.4 ±5.1 


7.8 ± 1.7 


3.7 ± 2.4 


N/A 


9.5 x 10- 


4 


3.2 ±0.9 ±1.0 


- 


B°-*D* a + TT- 






















2 




0.6 ±0.3 


< 0.14 


7.8 










D+^K*°K+ 


3 


2.8t?:S 


0.4 ±0.3 


0.3 ±0.2 


3.3 


3.9 x 10" 


2 


4 3+4-7 


< 12 


D+^K%K+ 







0.4 ±0.3 


< 0.14 


5.1 










all 


5 




1.2 ±0.4 


0.3 ±0.2 


N/A 


2.3 x 10" 


2 


1.911:3 ±0.5 


< 4.1 


B°^D7K+ 




















Df-*(/nT+ 


7 


5.8 ±2.6 


1.3 ±0.7 


1.1 ± 1.2 


13.0 


4.5 x 10" 


2 


2.4 ± 1.3 




D+^K*°K+ 


8 


7.3 ±2.9 


1.7 ±0.7 


< 0.7 


7.8 


1.9 x lO - 


3 


5.0 ±2.0 




D+^K° S K+ 


4 


3.7 ±2.0 


0.6 ±0.4 


1.3 ± 1.0 


9.2 


1.7 x 10" 


2 


2.5 ±2.1 




all 


19 


16.7 ±4.3 


3.5 ± 1.3 


2.7 ± 1.9 


N/A 


5.0 x 10" 


4 


3.2 ±1.0 ±1.0 




B°^Dt~K+ 



























0.8 ±0.6 


< 0.14 


5.3 










D+^K*°K+ 


1 




0.4 ±0.4 


< 0.14 


2.7 










D+^K° S K+ 


1 




0.4 ±0.4 


< 0.14 


4.3 










all 


2 




1.6 ±0.8 


< 0.14 


N/A 


0.48 






< 2.5 



fixed to the values obtained in a copious B ^D^*^ir Jt control sample. For the B°^D^ir~ and 
B°—>D~K + analyses, we obtain the threshold parameter £ from a fit to the data distributions of 
tties after loosening the M^ d and A.E requirements. In the case of B°—>D* + 7r~ and B°—>D*~K + , 
due to the low background level, we use simulated events to estimate £. 

No fit is performed to the B°^D*~ K + sample, due to the low number of events. Whenever 
there are enough events we perform a fit to each decay mode separately, as well as on the 
combination of all modes. The cross-contamination is estimated performing the same fit on the 
events in the data Mff nd sidebands (4a < \M^ d - 1968.6 MeV/c 2 | < 8a, where the resolution 
is a M cand = 5MeV/c 2 ). The number of observed events, the background expectations and the 
reconstruction efficiencies as estimated on simulated events are summarized in Table ffl. 

In the B°-kD+7t~ (B°^D~ K + ) mode the fit yields a Gaussian contribution of 21.4±5.1 (16.7± 
4.3) events and a combinatorial background of 7.8 ± 1.7 (3.5 ± 1.3) events. The cross-contamination 
is estimated to be 3.7± 2.4 (2.7± 1.9) events. The probability of the background to fluctuate to the 
observed number of events, taking into account both Poisson fluctuations and uncertainties in the 
background estimates, is 9.5 x 10 -4 (5.0 x 10 -4 ). For a Gaussian distribution this would correspond 
to 3.3cj (3.5a). Given the estimated reconstruction efficiencies we measure B(B°^>DfTr~) = (3.2 ± 
0.9) x 10~ 5 (B(B°^D-K + ) = (3.2 ±1.0) x 10~ 5 ), where the quoted error is statistical only. We also 
set the 90% C.L. limits B(B°^D* + ir-)< 4.1xl0~ 5 and B{B°^D*~K + )< 2.5xl0~ 5 . All results 
are preliminary. 

The systematic errors are dominated by the 25% relative uncertainty in B(Df — > 0tt + ). The 
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uncertainties on the knowledge of the background come from uncertainties in the £ parameter, 
for the combinatorial background, and from the limited number of events in the M^ d sidebands 
for the cross-contamination. They amount to 14%, 16%, 7% and 36% of the measured branching 
fractions in the B°^D+it-, B q ->D-R + , B ^D* + ir-&nd B°^D*~K + modes respectively. The 
rest of the systematic errors, which include the uncertainty on tracking, K® and charged kaons 
identification efficiencies range between 11% and 14% depending on the mode. 

In conclusion, we report a 3.3a signal for the b — > u transition £>°^.D+7r _ and a 3.5a signal for 
the B°^D~ K + decay, and we determine the preliminary results 



Since the dominant uncertainty comes from the knowledge of the Df branching fractions we also 
compute B(B°^D+tt-)xB(D+ -» 4>ir + ) = (1.13±0.33±0.21) x 1(T 6 and B(B°^DjK+) xB(Df -» 
= (1.16 ± 0.36 ± 0.24) x 10~ 6 . The search for B°-^D* + ir- and B°^D*-R + yields the 
preliminary 90% C.L. upper limits 
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